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Abstract 

Centrahty dependence of J/ip production in Au-Au collisions at y/sJiN = 200 GeV is studied 
using the Glauber model plus the kinetic formation model. Initial J/ip production and destruction 
by incoming nucleons are described by the Glauber model. J/'0-charm equilibration in the gluon 
plasma is described by the kinetic formation model. We explore the possibility of J/V' suppression 
suggested by the recent PHENIX data. We show that the J/V' yield monotonically decreases with 
increasing centrality in the kinetic formation model when the charm mass is smaller than a critical 
mass. The final J/ip yield is between the value after the Glauber suppression and the dynamical 
equilibrium value. This underscores the importance of both J ftp production from d-Au collisions, 
which is essential in determining the Glauber suppression, and high statistics Au-Au data, which 
can further constrain final state J/^-charm equilibration. 

PACS numbers: 25.75.-q,24.85.+p,14.40.Lb 



1 



I. INTRODUCTION 



Relativistic heavy ion colhsions have been used for the search of the Quark-Gluon Plasma, 
a deconfined nuclear matter that is believed to have existed at about one micro-second after 
the big-bang [l[0,0]. The existence and properties of the Quark-Gluon Plasma are reflected 
in the particle spectra produced in these collisions. One important particle that is sensitive 
to the Quark-Gluon Plasma production is the J/ip particle. As color screening in the Quark- 
Gluon Plasma can dissociate the J/ip particle, J/ip suppression relative to the yield from 
the superposition of nucleon-nucleon collisions can be used as a signal of the production of 
the Quark-Gluon Plasma 0|. 

In realistic collisions, the J / ip particle can also be destructed by incoming nucleons in the 
colliding nuclei. This is shown clearly by the proton-nucleus collision data where 
hot and dense nuclear matter is not expected to be produced. Collisions between the J/ij) 
particle and comovers produced in heavy ion collisions can further modify the J/ih yield 
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Hence detailed phenomenological models 



are necessary for the interpretation of J/ip production data. By comparing theoretical 
expectations and experimental results, an abnormal J/ip suppression in Pb-Pb collisions at 
the Super-Proton-Synchrotron (SPS) j22| has been identified which indicates the production 
of hot and dense nuclear matter at the SPS energies. 

At the Relativistic Heavy Ion Collider (RHIC) energies, there are on average more than 
one charm-anti-charm quark pair produced in each event. This opens the possibilit y o f 
producing the J/ ip particle in the deconfined phase via the reaction c-|-c— > J/'ip + g j23[ . 



Other possibi 
proposed 



ities involving J/ip equilibration at the hadronization time have also been 
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26l |. Many interesting studies have been made on the phenomenological 
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It was 



consequences of different production mechanisms 
demonstrated by the kinetic formation model that at collider energies, due to charm re- 
combination, J/ip enhancement may be the signal of Quark-Gluon Plasma production, in 
contrast to the traditional J/ip s uppression scenario. Recently, data on J/ip production at 
RHIC have been released j35|, |36[. No significant J/^ enhancement was observed. 

In this paper, we will explore the possibility of J/ip suppression in the framework of 
the kinetic formation model. In particular, we make the following variations of the original 
kinetic formation model. 
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'1) Recent pertubative Quantum Chromodynamics (QCD) studies of the charmonium 
spectrum give a charm mass m^^(m^^) = 1.243 GeV, different 

potential models |41| give = 1.5 GeV to = 1.8 GeV, preturbative QCD studies 
of charm production |42, |42| in hadron collisions give a charm mass from rrir = 1.3 
GeV to rUc = 1.5 GeV, a recent Nambu-Jona-Lasinio model calculation [2^ gives 
rric = 1.6 — 1.7 GeV. We will vary the charm quark mass between 1.3 GeV and 1.9 
GeV to study the effect of the threshold of J/ip production on J/ip-chaTm equilibration. 



44 



(2) J/^jj break-up cross section by a gluon was calculated in perturbative QCD 
The final state is a pair of charm and anti-charm quarks and the cross section decreases 
with the center-of-mass energy of the J/ ■j/'-gluon system. As the energy increases, the 
break-up is accompanied by additional parton radiation. We will model this effect by 
using a. J/ip break-up cross section that is a constant in energy and varies between 1 
mb and 6 mb and compare this with the dipole cross section results. 

(3) The original kinetic formation model uses an averaged uniform parton density in the 
transverse plane. This may underestimate charm-equilibration. In the following, we 
will use the Glauber model to study particle production and this leads to a better 
description of the evolution by local temperatures. 

(4) We make use of the recent discovery bv the PHENIX collaboration and various theo- 



retical approaches 



e recent discovery bv tne rrlJ 

Q,y,Q,Q,y,y,y,0 



that there is no strong charm energy 



loss in relativistic nuclear collisions. We also use the recent PHENIX J/ip data in pp 



collisions ^5| to normalize our calculations. 

In the next section, the Glauber model and the kinetic formation model will be introduced. 
We will use the Glauber model to describe the destruction of initial J/tp^s by incoming 
nucleons and use the kinetic formation model to describe J/ip equilibration in the Quark- 
Gluon Plasma. The results will be presented in section IIIII Then we will discuss the 
implications of our study. 



II. THE GLAUBER MODEL AND THE KINETIC FORMATION MODEL 



In nucleus-nucleus collisions, initial J/ip^s are produced by the collisions between incoming 
nucleons. Once a. J/ip is produced, other nucleons following the collision that produces the 

3 



i, y 




FIG. 1: Nucleus- nucleus collision transverse geometry. The reaction plane is determined by the 
X-axis that goes along the impact parameter vector b and the z-axis that goes along the beam 
direction. In this figure, the z-axis points out of the paper. The transverse position vector s starts 
at the center of the overlapping region. 

J/ ip will have opportunities to collide with it and turn it into a charm-anti-charm pair. The 
number of J/tp^s that are produced by collisions between incoming nucleons and survive 
further collisions with other nucleons can be written as 

NAB^j/^{h) = J dmNN^j/^tA{s + ^)tB{s- ^)Sa{s + ^)Sb{s- ^). (1) 

In the above formula, h is the impact parameter vector, and s is the transverse position. 
They are shown in Fig. [T] along with the transverse coordinate system. aMN^j/ip is the J/ip 
production cross section in nucleon-nucleon collisions. Ia and Ib are the nuclear thickness 
functions of nucleus A and nucleus B, while Sa and Sb are the average survival probabilities 
of a. J/tp after passing through nucleons in nucleus A and nucleus B. The normalization by 
the probability of a nucleus- nucleus collision at impact parameter b, <7AB{b), is not included 
as it only affects very peripheral events. 

We will focus on J/ip production in Au-Au collisions. To simplify the calculations, a 
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nucleus is taken to be a hard sphere with uniform nucleon density pn inside, and the radius 
is determined by i? = tqA^^^ {tq = 1.12 fm and A is the number of nucleons inside the 
nucleus). The number of survived J/ip^s can now be written as: 

where (Jj/^n is the J/ip break-up cross section due to the collision with an incoming 
nucleon and 2zo {2zq) is the longitudinal length of the row of nucleons inside nucleus 



A (B) at transverse position {x,y). More explicitly, zq = \J B? — {x + — y"^, and 
z'q = \J B? — {x — |)^ — y"^. The transverse position integration / ds becomes / dxdy over 
the region (x + |)^ + y"^ < and (x - + < ■ 

After the passing of the two incoming nuclei, the J/ip is immersed in the Quark-Gluon 
Plasma. As the initial gluon number is much larger than the initial quark number at RHIC 
energies js^ Is^, and in addition, light quark chemical equilibration time is much longer 
than the expansion time 5^ , the plasma will be approximated as a gluon plasma. As charm 



regeneration in the plasma is |57l . l58l |59| neglegible, charm quark number will be kept 
constant. There are indications that charm quark thermal equilibration may be possible 
[^, and we will assume local thermal equilibrium. 

J/ip equilibration in the plasma is achieved through the reaction J/ipg cc and its 
inverse reaction cc — > J/ipQ- Assume the rapidity distributions are fiat and focus on one 
fireball within one unit of rapidity, the J/tp number iVj/^ satisfies 

= \fN,p, - XnNj/^pg. (3) 

In the above, r is the proper time, Nc is the number of charm quarks per unit rapidity, pc 
is the density of anti-charm quarks and it equals the density of charm quarks, pc- Pg is the 
density of gluons. 

AiT' is the reduced J/ip formation rate through the reaction cc — > J/'ipO- It equals 
{o'cc^j/^gVrei,cc) , whcrc (■ ■ ■) dcuotcs average over charm and anti-charm quarks, and Vrei,cc 
is the relative speed between a charm quark and an anti-charm quark. When the evolution 
is in thermal equilibrium, i.e., the evolution is controlled by one temperature T, 

= or^ I ^^o / ^A^^ ' ^al)z' K,{z)a,,^j/^g{s = z'T^). (4) 
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We use ac to stand for the mass of charm quark divided by the temperature, i.e., ttc = tric/T. 
zq = max(Q;c + ac,aj/^ + ag). Kn{x) is the nth order modified Bessel function, s is the 
center of mass energy squared. 

Similarly, A/?, the reduced J/ip destruction rate via the reaction J/ipg cc can be 
written as 

1 f°° 

= 7^/^ V / ^^(^^ - a]/^fKi{z)aj/^g^ccis = z^T^). (5) 

The detailed balance relation, 

acc^j/^g{s) = — ——-^aj/^g^cc[s), [b) 

(^cc yf^cc) 

relates the two reduced rates. In Eq. dj/^g = 48 is the number of degrees of freedom 
for the J/ipg system, and d^c = 36 is the number of degrees of freedom of the cc system. 
kj/tpg and kcc are the center-of-mass momenta of the J/ipg and cc systems. After using the 
detailed balance relation, we can rewrite Xp as. 



= -d- i2alK,ia,)y L " "^^^^ K,{z)aj/,g^,,{s = z T ) 

dj/^g 2«j/^^2(aj/v.) 

In Eq. 1^, Nc = Ncfi — Nj/^ where Ncfi is the number of charni quarks per unit rapidity 
at the initial time tq. = Nc/V{t). We use the Bjorken model 61| for the evolution and 
the volume V{t) = Vqt/to. The gluon density pg = ^C{3)T^ ^ l.gST^ = 1.95^. 

III. J/ij PRODUCTION AT ^5]^ = 200 GEV 
A. Initial conditions 

The initial gluon number, charm-anti-charm quark number, and J/ip number for central 
(6 = 0) Au-Au collisions will be specified. For non-zero impact parameters, the total number 
of gluons and the total number of charm- ant i- charm quarks follow binary collision scaling, 
and the J/ip number follows Eq. (jll2|) . The transverse density per unit rapidity of gluons 
(charm-anti-charm pairs) at impact parameter b and position s is proportional to ^^^(^'-l- 
- I) and is given by 



where T^uAuib = 0) is the nuclear thickness function at zero impact parameter. It equals 
JdstAu{s + l)tAu{s- I) = AplJds{R^ - s^) = with A = 197 for Au nucleus. The 

initial temperatures are local and calculated from the initial gluon density distribution. The 
J/ip transverse density per unit rapidity needs additional correction factor SAuis + ^)SAu{s — 
|). It can be written as 

^^(6, s) = ^PNZ^^'o j.^^2{b = Q)^^dy^ ^^ ^ 0)F{2pNaj/^r,zo)F{2pNaj/^Mz'o), (9) 

where function F(x) = Iz^EEizii, The initial time is taken to be Tq = 0.2 fm/c. In central 
collisions, the initial gluon rapidity density is = 300 and charm rapidity density is 
= 2.5. The y a re in line with perturbative Quantum Chromodynamics (QCD) based 



63|. The initial J/ip rapidity density before correcting for the Glauber 



dy 

calculations _ 

suppression is ^^^^ = 0.033. This value comes from the superposition of J/ip production 
from nucleon-nucleon collisions at ^/SNN = 200 GeV using the production cross section 
measured by the PHENIX collaboration |3i| and the nuclear thickness function at zero 
impact parameter. The J/ip mass is fixed at rrij/^ = 3.1 GeV. The charm quark mass 
will be varied between 1.3 GeV and 1.9 GeV. The freeze-out temperature Tf is taken to be 
0.15 GeV. 



B. Cross sections 
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nucleon dissociation cross section aj/^N will be varied among 0, 4.4 mb, and 



The JM-nucleon dissoc 
1 mb Q, n, 64, 65, Q 



651 166|. The J/ipg cc cross section is taken to be a constant. We 
will vary (Tj/^g^cc to study J / ip-chaxm. equilibration. A comparison with the results from a 
dipole cross section will be made to study the effect of the energy dependence of the cross 
section. 

Fig. (0) shows the J/ip production and destruction cross sections as functions of the 
center-of-mass energy squared s. In Fig.'s (I2t)-(I2t)5 crj/'tpg^cc is kept at a constant value 
of 3 mb, which is about the J/ip size in the vacuum. The inverse cross section, (Jcc^j/ipg, 
depends on the charm quark mass. When the charm quark mass ttIc = 1.5 GeV, which 
is smaller than the J/ip production threshold, the inverse cross section starts at zero and 
increases with s. It crosses the destruction cross section at s ~ 40 GeV^. However, when 
nic = 1.7 GeV, which is larger than the J/ip production threshold, there is a peak at the 
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FIG. 2: The J/^l^ production and destruction cross sections in a gluon plasma as functions of the 
center-of-mass energy squared s. In (a)-(c), crj/^g = 3 mb. In (d), crj/^pg is specified in Eq. (^Oil- 
threshold. When s increases above the threshold, acc^j/^g first decreases to below aj/^g^cc- 
It then increases and crosses Ojj^g^cz at s ~ 30 GeV. The width of the peak at the threshold 
depends on the charm quark mass. The larger the mass, the wider the peak. This can be 
seen by comparing a^c^jj^g for rric = 1.7 GeV in Fig. ^p) with a^c^ji^g for rric = 1.87 GeV 
in Fig. (j2fc)- In addition, when rric increases, the crossing position decreases. Eventually, 
acc^j/^g is always larger than ojj^g^cc as shown in Fig . (Eh) . We also plot the Jli)g — >■ cc 
cross section based on a perturbative QCD calculation |44l. l45l| . 



and its inverse cross section in Fig. (I21i)- In the above equation, yU = mc/2 is the 
reduced charm quark mass in the charm-anti-charm center-of-mass system, k is the gluon 
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FIG. 3: The reduced reaction rates, {(Tjj^g^ccV) and {(Jcc^j/ipgv) , as functions of the temperature 
of the plasma. 

momentum, and cq = — mj/^ is the binding energy of J /ip. For mc equals the D meson 
in vacuum mass 1.87 GeV, the J /if) destruction cross section has a finite peak near the 
threshold. The peak of the cross section is about 2 mb. The production has an infinite peak 
at the threshold, and it is always larger than the destruction cross section. 

The reduced J/ip production rate {pcc^ji^gv) and destruction rate {p ji^g^ccV) for differ- 
ent cross section choices are plotted in Fig. 0. For smaller than the J/V^ production 
threshold, when aj/^g^cc is a constant, the reduce J/ip destruction rate is always a constant 
and it equals o-jj^g^ccC- In Fig.'s (jHtL)-©), (yj/^g-^cc = 0.3 fm^. For nic = 1.5 GeV, which is 
smaller than the J/ip production threshold, the destruction rate is a constant and the pro- 
duction rate increases with the temperature. The destruction rate is always larger than the 
production rate. As increases, the production rate increases. When is larger than the 
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J/if) production threshold, the destruction rate decreases with increasing iric. For rric = 1.7 
GeV, {(Xji^g^ccv) starts from a value smaller than aj/^g^ccC- It increases with temperature 
to approach aj/^g^ccC = 0.3 fm^c. {<Jcc^j/^gv) has the same shape as the nic = 1.5 GeV 
case and increases with the temperature of the plasma. It is smaller than the destruction 
rate up to T = 0.7 GeV. When = 1.87 GeV, which is larger than the production thresh- 
old, the production rate is always larger than the destruction rate. The production rate for 
TTLc = 1.87 GeV decreases as the temperature of the plasma increases, reflecting the decrease 
in the production cross section away from the threshold. We also plot the rates for the case 
with the dipole J/^ipg — cc cross section determined by Eq. (|Tn|l in Fig. (jHl). ((Xcc^j/^gf) is 
still larger than {<Jj/^g^cc'v)- However, because the peak value of the dipole cross section is 
smaller than 0.3 fm^, also because the cross section is exclusive and decreases with s, the 
rates are smaller than the constant cross section case shown in Fig. (jSt). 



C. Centrality and normalization 

We will use the number of participants to characterize the centrality. The number of 
participants is calculated in the Glauber model with hard sphere geometry according to 



ip) = J ds{tA{s + -)(1 - exp{-aNNtB{s- -))) 
b b 

+tB{s- -)(1 - exp{-(TNNtA{s + -)))) 



= 8pN dx dy{zoil - exp{-2(TNNpNZo)) 

+Zq{1 - exp{-2aNNPNZo))), (11) 

In the above, aNN is the nucleon-nucleon inelastic cross section. We will use cnn = 40 mb. 

The J/ Ip central rapidity density is normalized to the number of binary nucleon-nucleon 
collisions. We will plot B '^^J^''' \ y=o/Ncoii, where B = 0.06 is the J/ip to e~^e~ or decay 
branching ratio, the number of binary nucleon-nucleon collisions is given by 

b. _ 6, 



f b b 

Ncou{b) = (Jnn J dstA{s + -)tB{s- -] 



= 4:p^aNN / dx / dyzoZQ. (12) 

Jo Jo 

The number of participant nucleons and the number of binary nucleon-nucleon collisions 
are plotted in Fig. (jH). 
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FIG. 4: The number of participant nucleons Npart and the number of binary nucleon-nucleon 
colUsions Ncoii as functions of the impact parameter. 

D. Evolution of the gluon plasma 

The evolution of the gluon plasma is described by the Bjorken model. We will use the 
latest freeze-out time tj to characterize the lifetime of the plasma. Fig. (jH)) shows that 
as the number of participant nucleons increases, the lifetime of the plasma increases. The 
typical lifetime for a plasma with an initial gluon rapidity density of 300 and a freeze-out 
temperature of 0.15 GeV is about 5 fm/c. Increasing the gluon rapidity density to 500 will 
bring the lifetime to about 8.5 fm/c, while increasing the freeze-out temperature to 0.17 
GeV will decrease the lifetime to about 3.5 fm/c. 

The temperature in the central transverse cell at the initial time as a function of the 
number of participating nucleons is shown in Fig. ©. It increases very rapidly when the 
number of participating nucleons goes from to 50 and changes only slowly when the number 
of participating nucleons goes above 100. The typical initial temperature in the central cell 
is about 500 MeV. 
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FIG. 5: Freeze-out time Tt as a function of the number of participant nucleons Npart- 
E. J production results 

J/il) enhancement was predicted by Thews et. al using the kinetic formation modeL The 
PHENIX data excluded strong enhancement of the J/ip particle. We are going to explore 
the possibility of J /ip suppression from the kinetic formation model. About 10 cc pairs are 
produced on average in each Au-Au event. We assume no strong initial state energy loss of 
the produced cc pairs. The rapidity range of the cc pairs is taken to be Ay = 4. Now we 
vary the charm quark mass mc, and observe the effect of changing the dynamical equilibrium 
of the J / ip particle through the reactions cc J / ipg and J/ ipg — cc. 

The destruction of J/ip or its precursor (the correlated cc pair that produces a J/ip 
when there are no final state interactions) by collisions with incoming nucleons depends on 
the cross section crj/^i^. When the colliding nuclei are highly Lorentz contracted at high 
enough energies, the nuclei pass through each other much earlier than the formation time 
of a J/ip- The J/ip particle will not be strongly suppressed by the incoming nucleons. We 
will use (Tj/^jy = to simulate this case. The centrality dependence of J/ip yields are 
shown in Fig. When mc=1.3, 1.5, 1.7 GeV, the final J/ip number is smaller than the 
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FIG. 6: Initial temperature of the center cell as a function of the number of participant nucleons 

superposition of nucleon-nucleon collisions. The larger the cross section, the smaller the 
final J/ip number. When the J/ip destruction cross section crj/^g goes up over about 5 
mb, the curves converge. This gives the approximate dynamical equilibrium J/ip value. 
The equilibrium value appears to be independent of centrality for mid-central and central 
Au-Au collisions. It increases with the charm quark mass rric- Also notice that rric — 1.7 
GeV is larger than the J/t/j production threshold. However, the initial J/t/j is higher than 
the equilibrium value, and the final J/ip number is smaller than the initial. For rUc = 1.9 
GeV, The equilibrium value is larger than the initial J/i/j number and the final J/ip number 
increases with increasing J/ip gluon cross section. In the J/ip suppression case, the final J/ip 
number decreases with increasing centrality while in the J/ ip enhancement case, it increases 
with increasing centrality. 

The adaptive step size Runge-Kutta method is used in solving the evolution equation. 
The transverse plane is divided into 20x40 cells (20 divisions in the x direction from — (i?— |) 
to i? — |, and 40 divisions in the y direction from —\JR'^ — (|)^ to -^^^"-^"(1)^) . Results are 
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FIG. 7: Centrality dependence of J/ip central rapidity density per binary nucleon-nucleon collision 
in Au-Au collisions when the suppression of J/ip by incoming nucleons is negligible {ctj/^n = 0)- 
Solid lines are for 20 x 40 cells and dashed lines are for 10 x 20 cells. 

shown in Fig. ((7j) as solid lines. Dashed lines are for the 10 x 20 cell division case. They 
agree well with the 20 x 40 cell division results. 

At RHIC energies, the number of J/ip is much smaller than the number of charm-anti- 
charm quark pairs. One can approximate the number of charm quarks at time r by the 
initial number of charm quarks, and solve Eq. © to obtain the solution 

f"^! \F{j)dT 



Nj/A'^f) = e(r/,i-o) 



TO V{T)e{T,To) 



(13) 



Here, 



e(r. To) = exp |- j^^^ AB(r')p,(r')f/r'| . (14) 

The approximation solution gives no observable difference when the J/ip production rate is 
small, for example, when rUc = 1.3 GeV. When the production rate is large, for example. 
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FIG. 8: J ftp number evolution for central (6 = 0) Au-Au collisions at y^s/vTV = 200 GeV. 

when rric = 1.9 GeV and Cj/^g = 6 mb, it can overestimate the J/ip yield (relative to the 
exact solution) by about 4%. 

We show the time evolution of J/ip number in central {b = 0) collisions in Fig. (jHJ. 
When rric =1.3, 1.5, 1.7 GeV, the J/ip number changes drastically before r =1 fm/c. When 
mc=1.7 GeV and crj/^g =5, 6 mb, the J/ip value reaches the lowest point at about 1 fm/c, 
then it increases slightly. This is because of the production of J/ip in central cells. For 
rric = 1-9 GeV, the J/ip number changes over the entire evolution. The equilibration time is 
larger than the system lifetime and results with different cj/^g values do not converge and 
dynamical equilibrium of J/ if) has not been reached over the evolution. 

To check the sensitivity on parameters, we change the gluon rapidity density from 300 
to 500. Results are shown in Fig. ©. As the initial gluon number increases, the initial 
temperature increases, more J/ip's will be destructed during the evolution. This leads to an 
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FIG. 9: Centrality dependence of J/ij: central rapidity density per binary nucleon-nucleon collision 
when (Tji^jq = mb and dNg/dy{h = 0) = 500. 

overall decrease of the J/ip yield for all charm mass values. If we keep the gluon rapidity 
density to be 300 and increase the freeze-out temperature from 0.15 GeV to 0.17 GeV, no 
significant changes show up for ttIc =1.3,1.5, 1.7 GeV as seen from Fig. (fTIUl . This indicates 
that the equilibration time is much longer than the time that goes from T = 0.17 GeV 
to T = 0.15 GeV. For the rric =1.9 GeV case, especially for large <7j/^g values, the final 
J/ip values decrease significantly. For example, when aj/^g = 6 mb, the final J/ip value 
for central collisions decreases from 0.31 to about 0.22. This shows that for the rric = 1.9 
GeV case, the final J/tp value is still far from equilibrium and the equilibration time is not 
negligible compared with the time that goes from T = 0.17 GeV to T = 0.15 GeV. 

To see the difference between the constant aj/^g^cc case and the dipole form aj/^g^cc 
specified in Eq. f|TU|) . we show the centrality dependence of J/iIj production in Fig. (fTT|l. 
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FIG. 10: Centrality dependence of J/ip central rapidity density per binary nucleon-nucleon collision 
when (Tji^jq = mb and Tj = 0.17 GeV. 

The J /if) centrality dependence in the dipole cross section case follows closely the constant 
cross section case. It is between the two curves with <Jj/^g^cc = 1-0 nib and ajj^g^cc = 2.0 
mb. The peak value of the dipole cross section is about 2 mb. The J /ip yield is smaller than 
the constant 2 mb destruction cross section case because the dipole cross section does not 
have destruction by high energy gluons with additional radiation involved. The enhancement 
relative to the superposition of nucleon-nucleon production when aj/^g^cc = 2.0 mb is about 
1.2 times larger than the enhancement when the dipole cross section is used, independent 
of centrality. In other words, the dipole cross section case can be described by the constant 
cross section case with an effective K factor. We notice that the dipole cross section and the 
constant cross section have very different dependences on the center-of-mass energy (Fig.|2I). 
The above comparison indicates that the centrality dependence of J/ip production reflects 
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FIG. 11: Centrality dependence of J/ip central rapidity density per binary nucleon-nucleon collision 
when for rric =1.87 GeV. 

the overaU production rate and can not differentiate between detailed cross section shapes. 

Now we add in the Glauber suppression of J/ip- We will first look at the (Tj/^n = '^■'^ nib 
case. For =1.3, 1.5, and 1.7 GeV, the final J/ip values are still smaller than the initial 
values as shown in Fig. ()12j). Because of the destruction by incoming nucleons, the final 
J/ip values and the initial are closer. For ttIc = 1.9 GeV, the final J/tp values first decrease 
with centrality and then increase with centrality. This centrality dependence differs from 
the monotonic centrality dependence when rric = 1.3, 1.5, 1.7 GeV. Because of suppression 
by the incoming nucleons, when crj/^g is small, the final J/ip values are no longer larger than 
the superposition of J/tp from nucleon-nucleon collisions. If aj/^N is further increased to 7.1 
mb, the initial J/ip number after collisions with incoming nucleons is even smaller as shown 
in Fig. p3p . The range between the initial and final J/ip values decreases for the rric =1.3, 
1.5 GeV cases. For the = 1.7 GeV case, it goes below the equilibrium value and the final 
state J/ip equilibration increases the J/ip yield. When rUc = 1.9 GeV, the range between 
initial and equilibrium increases. 
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FIG. 12: Centrality dependence of J/ if) central rapidity density per binary nucleon-nucleon collision 
when (Tjj^j^ = 4.4 mb. 

IV. SUMMARY AND CONCLUSIONS 

The J /ip production in Au-Au collisions at y/sJm = 200 GeV is studied in the framework 
of the kinetic formation model. It is shown that when the charm quark mass rric is smaller 
than some critical value, the J/ip yield monotonically decreases with increasing centrality. 
As (Tjj^g-^cc increases, the final J/%1) yield moves away from the value right after the Glauber 
suppression by the incoming nucleons and approaches the dynamical equilibrium value. The 
J/il) production data from d-Au collisions at the Relativistic Heavy Ion Collider will help 
to fix the Glauber suppression and the higher statistics Au-Au data will further help to 
understand J / ^-charm equilibration. 

The current model does not have gluon nuclear shadowing. Further inclusion of shadowing 
will decrease the number of charm quarks produced and lower the number of J/'^'s produced. 
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FIG. 13: Centrality dependence of J/ip central rapidity density per binary nucleon-nucleon collision 



when a 



J/ipN 



7.1 mb. 



Color screening j4 , 



67 



68, 



69| will lead to more dissociation of J/ and more suppression of 



final J/ ip yield. The inclusion of feeddown will lead to the same effect as higher charmonium 
resonances are easier to be dissociated. Transverse expansion will lead to earlier freeze-out 
and should not change the qualitative feature of J / tp suppression. J / ^/^-charm equilibration 
in the hadron phase is beyond the scope of this paper. However, as the equilibrium value 
in the hadron phase is smaller than the superposition of production from nucleon-nucleon 
collisions, equilibration in the hadron phase will also lead to J/ip suppression in Au-Au 
collisions. 



20 



Acknowledgments 



We thank B.A. Li and Z. Lin for helpful discussions. We also thank B. Hemphill and T. 
Pranks for participation in the early stage of the study. This work is supported by the U.S. 
National Science Foundation under Grant No. PHY-0140046, and the Arkansas Science and 
Technology Authority under Grant No. Ol-B-20. 



[1] B. Miiller, Rep. Prog. Phys. 58, 611 (1995). 

[2] S.A. Bass, M. Gyulassy, H. Stocker, and W. Greiner, J. Phys. G 25, Rl (1999). 

[3] Proceedings of 16th International Conference on Ultra-Relativistic Nucleus-Nucleus Collisions 

(Quark Matter 2002), Nucl. Phys. A715. 
[4] T. Matsui and H. Satz, Phys. Lett. B 178, 416 (1986). 
[5] E772 Collaboration (D.M. Aide et al), Phys. Rev. Lett. 66, 133 (1991). 
[6] E866/NuSea Collaboration (M. Leitch et al), Phys. Rev. Lett. 84, 3256 (2000). 
[7] NA 50 Collaboration (M.C. Abrcu et al), Phys. Lett. B 410, 337 (1997). 
[8] E.V. Shuryak, Sov. J. Nucl. Phys. 28, 408 (1979). 

[9] X.M. Xu, D. Kharzeev, H. Satz, and X.N. Wang, Phys. Rev. C 53, 3051 (1996). 
[10] S. Gavin and R. Vogt, Phys. Rev. Lett. 78, 1006 (1997). 

[11] D. Kharzeev, C. Lounrenco, M. Nardi, and H. Satz, Zeit. Phys. C 74, 307 (1997). 

[12] W. Gassing and CM. Ko, Phys. Lett. B 396, 39 (1997). 

[13] C. Gale, S. Jeon, and J. Kapusta, Phys. Lett. B 459, 455 (1999). 

[14] J. Geiss, C. Greiner, E.L. Bratkovskaya, W. Gassing, and U. Mosel, Phys. Lett. B 447, 31 
(1999). 

[15] D.E. Kahana, S.H. Kahana, Phys. Rev. C 59, 1651 (1999). 

[16] B.H. Sa, A. Tai, H. Wang, and F.H. Liu, Phys. Rev. C 59, 2728 (1999). 

[17] C. Spieles, R. Vogt, L. Gerland, S.A. Bass, M. Bleicher, H. Stocker, and W. Greiner, Phys. 

Rev. C 60, 054901 (1999). 
[18] R. Vogt, Phys. Rep. 310, 197 (1999). 

[19] J.P. Blaizot, P. M. Dinh, J.Y. Ollitrault, Phys. Rev. Lett. 85, 4012 (2000). 

[20] A. CapeUa, E.G. Ferreiro, and A.B. Kaidalov, Phys. Rev. Lett. 85, 2080 (2000). 



21 



[21] C. Gale, S. Jeon, J. Kapusta, Phys. Rev. C 63, 024901 (2001). 

[22] NA50 Collaboration (M.C. Abreu et al), Phys. Lett. B 477, 28 (2000). 

[23] R.L. Thews, M. Schroedter, and J. Rafelski, Phys. Rev. C 63, 054905 (2001). 

[24] P. Braun-Munzinger and J. Stachel, Phys. Lett. B 490, 196 (2000). 

[25] M.I. Gorentein, A.P. Kostyuk, H. Stocker, and W. Greiner, Phys. Lett. B 509 277 (2001). 

[26] L. Grandchamp and R. Rapp, Phys. Lett. B 523, 60 (2001). 

[27] B. Zhang, CM. Ko, B.A. Li, Z. Lin, B.H. Sa, Phys. Rev. C 62, 054905 (2000). 

[28] M.I. Gorentein, A.P. Kostyuk, H. Stocker, and W. Greiner, Phys. Lett. B 524 265 (2002). 

[29] L. Grandchamp and R. Rapp, Nucl. Phys. A 709, 415 (2002). 

[30] B. Zhang, CM. Ko, B.A. Li, Z.W. Lin, and S. Pal, Phys. Rev. C 65, 054909 (2002). 

[31] B. Zhang, B.A. Li, A.T. Sustich, and C. Teal, Phys. Lett. B 546, 63 (2002). 

[32] A. Andronic, P. Braun-Munzinger, K. Redlich, J. Stachel, Phys. Lett. B 571, 36 (2003). 

[33] E.L. Bratkovskaya, W. Gassing, H. Stocker, Phys. Rev. C 67, 054905 (2003). 

[34] L. Grandchamp, R. Rapp, G.E. Brown, e-Print archive: hep-ph/0306077 

[35] PHENIX Collaboration (S.S. Adler et al), Phys. Rev. Lett. 92, 051802 (2004). 

[36] PHENIX Collaboration (S.S. Adler et al), Phys. Rev. C 69, 014901 (2004). 

[37] N. Brambilla, Y. Sumino, and A. Vairo, Phys. Lett. B 513, 381 (2001). 

[38] M. Eidemiiller and M. Jamin, Phys. Lett. B 498, 203 (2001). 

[39] S. Recksiegel and Y. Sumino, Phys. Rev. D 67, 014004 (2003). 

[40] S. Rechsiegel and Y. Sumino, Phys. Lett. B 578, 369 (2004). 

[41] E.J. Eichten and C. Quigg, Phys. Rev. D 49, 5845 (1994). 

[42] R.V. Gavai, S. Gupta, P.L. McGaughey, E. Quack, P.V. Ruuskanen, R. Vogt, and X.N. Wang, 

Int. J. Mod. Phys. A 10, 2999 (1995). 

[43] G. Altarelh, M. Diemoz, G. Martinelh, and P. Nason, Nucl. Phys. B308, 724 (1988). 

[44] M.E. Peskin, Nucl. Phys. B156, 365 (1979). 

[45] G. Bhanot and M.E. Peskin, Nucl. Phys. B156, 391 (1979). 

[46] PHENIX Collaboration (K. Adcox et al), Phys. Rev. Lett. 88, 192303 (2002). 

[47] Y.L. Dokshitzer and D.E. Kharzeev, Phys. Lett. B 519, 199 (2001). 

[48] S. Batsouh, S. Kelly, M. Gyulassy, and J.L. Nagle, Phys. Lett. B 577, 26 (2003). 

[49] M. Djordjevic and M. Gyulassy, Phys. Lett. B 560, 37 (2003). 

[50] M. Djordjevic and M. Gyulassy, Phys. Rev. C 68, 034914 (2003). 

22 



[51] B.W. Zhang, E. Wang, and X.N. Wang, e-Print archive: |nucl-th/ 0309040| 

[52] M. Djordjevic and M. Gyulassy, Nucl. Phys. A733, 265 (2004). 

[53] D. Kharzeev and K. Tuchin, Nucl. Phys. A735, 248 (2004). 

[54] A.D. Martm, R.G. Roberts, W.J. Stirling, and R.S. Thorne, Eur. Phys. J. C 23, 73 (2002). 

[55] K.J. Eskola, Nucl. Phys. A590, 383c (1995). 

[56] X.N. Wang, Phys. Rept. 280, 287 (1997). 

[57] P. Levai, B. Miiller, and X.N. Wang, Phys. Rev. C 51, 3326 (1995). 

[58] Z. Lin and M. Gyulassy, Phys. Rev. c 51, 2177 (1995). 

[59] Z.W. Lin and M. Gyulassy, Nucl. Phys. A590 495c (1995). 

[60] B. Svetitsky, Phys. Rev. D 37, 2484 (1988). 

[61] J.D. Bjorken, Phys. Rev. D 27, 140 (1983). 

[62] X.N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 (1991). 

[63] M. Gyulassy and X.N. Wang, Comput. Phys. Commun. 83, 307 (1994). 

[64] J.L. Nagle and M.J. Bennett, Phys. Lett. B 465, 21 (1999). 

[65] M.J. Bennett, J.L. Nagle, Phys. Rev. C 59, 2713 (1999). 

[66] NA50 Collaboration (L. Ramello et al), Nucl. Phys. A715, 243 (2003). 

[67] S. Datta, F. Karsch, P. Petreczky, and I. Wetzorke, Nucl. Phys. Proc. Suppl. 119, 487 (2003). 

[68] T. Umeda, H. Matsufuru, O. Miyamura, and K. Nomura, Nucl. Phys. A721, 922 (2003). 

[69] M. Asakawa and T. Hatsuda, Phys. Rev. Lett. 92, 012001 (2004). 



23 



